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The coupled line equations for two modes traveling in the same direction 
are considered. The covariances of the mode transfer functions are cal- 
culated in the case of random coupling. Exact results are obtained when 
the coupling is a function of a finite state Markov chain, and also when 
the coupling is white noise. Perturbation results are obtained in the case 
of weak, zero mean, wide sense stationary coupling. It is also shown that 
perturbation results are valid in the case of strong coupling, if the cor- 
relation length is short. 

I. INTRODUCTION 

We consider the coupled line equations for two modes traveling in 
the forward direction, 1 



^ + r / (2) = JC(Z)J>(2), 

^ + r,/,(z) = jc(z)i (z), 



(1) 



subject to the initial conditions 

/o(0) = i , 7,(0) = i t . (2) 

The time dependence exp (2irjft) has been suppressed. Here c(z) is 
the real coupling coefficient, which we take to be a random function 
of z, and the loss and phase constants are 

To = «o + ?/So , I\ = «! + jpi . (3) 

These equations provide an approximate description of a variety of 
physical systems 2 ' 3,4 such as optical fibers 5,6 and metal waveguides. 7,8 
Typically, with the choice i a = 1, i, = 0, 7 (z) represents a desired 
mode launched at z = 0, and 7i(z) an undesired spurious mode. 
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Rowe and Young 1,9 have studied these equations with c(z) taken 
to be white noise. When c{z) is white noise, equations (1) are symbolic, 
and Rowe and Young interpret them by a limiting process, and with 
the aid of a matrix technique obtain exact expressions for the means 
and covariance functions of the solutions. With the aid of these solutions 
they studied pulse distortion caused by the coupling, and provide 
additional verification of Personick's 10 interesting result that in some 
cases increasing the random coupling can improve the impulse response. 

This article has three main purposes. The first is to show that there 
is another class of stochastic coupling coefficients for which the means 
and covariances of the solutions of (1) can be calculated exactly. These 
involve Markov processes which can assume only a finite number of 
values. We set up the equations satisfied by the covariance functions 
and means for a general coupling of this type in Section II, and solve 
them in a special case in Section III. 

The second purpose of this paper is to show, with the use of Ito 
integrals and the Ito calculus, that the interpretation of equations (1) 
by Rowe and Young when c(z) is white noise is consistent with the 
interpretation due to Stratonovich." This is the subject matter of 
Section IV. 

The third purpose is to apply a perturbation technique developed 
by Papanicolaou and Keller 12 to the coupled line equations. It is shown 
that this method can be applied not only to the case of weak coupling, 
but also to the case of strong coupling if the correlation length is short. 
This is the subject of Section V. 

In the remainder of Section I, we define the quantities of interest, 
introduce notation, and discuss some of the previously mentioned 
points at greater length. 

The differential loss and phase constants are given by 

Ar = To - T, = Aa + ;'A/3. (4) 

We write the coupling coefficient in the form 

c(z) = CN(z), (5) 

where N(z) is dimensionless. Following Rowe and Young 1 '" we assume 
that Ac* is independent of the frequency /, but that A/3 and C are odd 
functions of /. We further assume that C is an odd function of A/3. 
Also, let. 

g (z, A/3) - e r °*h(z), < h {z, A/3) = e v °'I>(z). (6) 

Note that our g n (z, A/3) corresponds to Rowe and Young's (!,, (z) and 
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our g L (ft, A/3) to their G,(«). Then, from (1), (4), and (5), 

(7) 



*S = jC(W)N(z) gi (z, A/3) 



^ = Ar ?1 ( 2 , A/3) + jC(W)N(z)g (z, A/3), 

with initial conditions 

<7o(0, A/3) = io , 0.(0, A/3) = i, . (8) 

Now define the correlation functions 

Rkiiz) = iaifyh W + ■)**<* A8)>, (fc, I = 0, 1), (9) 

where * denotes complex conjugate and ( ) denotes stochastic average. 
Rowe and Young 1,9 have calculated R kt {z) exactly for the case of white 
noise coupling, with spectral density D , 

(N(z)N(fi) = D 8(z - J), (10) 

in the case of initial values i" = 1 and i t = 0, corresponding to the 
signal and spurious modes respectively. They have also calculated 
the average of the squared envelope of the impulse response, which 
involves n double integral with respect to A/3 -and a of R 00 (z). The 
calculation was exact for frequency-independent coupling' (i.e., C = 
c sgn A/3), and approximate for moderate fractional bandwidths for 
frequency-dependent coupling. 9 

In this paper we will be concerned with the calculation of R kc {z) 
only, and will not consider the time domain statistics. We note that 
the average powers in the two modes are given by 

I\(z) = e- ta "R 00 (z) |,-o. l\(z) = e- 2a ° ! R u (z) |„. . (11) 

We first consider the case 

N(z) = F(M(z), z), (12) 

where M{z) is a finite state Markov chain 13 which has, in general, 
a nonstationary transition mechanism. Using the results of a paper 
by Morrison, 14 we obtain a system of ordinary differential equations, 
with prescribed initial conditions, for calculating R ki (z). If the process 
]\I{z) has a stationary transition mechanism, and F is a function of 
M(z) alone, these equations have constant coefficients. In the particular 
case in which N(z) = T(z), where T(z) is the random telegraph process, 15 
we obtain the Laplace transforms of R k( explicitly. 
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Next we investigate the case of white noise coupling by means of 
the Ito calculus. 16 This approach differs from that used by Rowe and 
Young, 1 ' 9 who consider the line as the limit of discrete sections, of 
vanishing length, of uncoupled ideal lines, with discrete mode converters 
at the end of each section. The quantities R ke (z) are calculated exactly 
for white noise coupling, in agreement with the results of Rowe and 
Young. 1,9 

Finally, we turn to the asymptotic calculation, for long lines, of 
the quantities R kl (z), in the case of weak coupling, weak attenuation, 
and narrow fractional bandwidth. The coupling is assumed to have 
zero mean and to be wide sense stationary, so that c(z) is given by 
(5), where 

(N(z)) = 0, <iV(z)AT(f)> = p(z - f). (13) 

The asymptotic equations for R kt (z) are determined, first for the 
nonresonance case, and then for the resonance case, in which the 
differential phase constant is small. There is an alternate interpretation 
of the resonance case, for which the correlation length is short. It turns 
out that for white noise coupling (which has zero correlation length), 
corresponding to (10), the asymptotic equations for R Q0 (z) and R n (z) 
are exact, in both the nonresonance and resonance cases, and those for 
R {) i(z) and R l0 (z) are exact in the resonance case. 

II. MARKOV CHAIN COUPLING 

We begin by writing down the equations satisfied by the quantities 

r kl {z) = g k (z, A/3 + cr)r/*(z, A/3), (k, t = 0, 1). (14) 

Let 

Co = C(A/3), C„ = C(A/3 + «r). (15) 

Then, from (4) and (7), for sufficiently smooth N(z), it follows that 
^ = jN(z)(C,r u , - (V 0I ), (16) 

^ = jN(z)(C,r n - 0„„) + Ar*r,„ , (17) 

^ = jN(z)(Cs oa - Cnru) + (Ar + ja)r l0 , (18) 

dz 

^i = jN(z)(C a r in - C o r 10 ) + (Ar* + Ar + j<r)r u . (19) 
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These equations do not hold when N(z) is white noise, as will be dis- 
cussed later. From (8) the initial conditions are 

r*,(0) = i k i* t , (k, t = 0, 1). (20) 

We now consider the case of coupling corresponding to (12), and 
first state some properties of the finite state Markov chain. 13 The 
sample functions M{z) are defined on the half line ^ z < °o , have 
right-continuous paths, and can take on only a finite number N of 
distinct values a p (p = 1, ■ • ■ , .V). An initial probability distribution 
is given: 

X„ = Prob {;l/(0) = a p ], (p = 1, •■■ , N), (21) 

where X p > and 

ZX„=1. (22) 

P =i 

The transition probabilities are defined, for ^ x ^ y, by 

P„(x, y) = Prob \M{y) = a, \ M(x) = a p \, 

(p, q = 1, • • • , N). (23) 

We consider only those processes which can be defined by means 
of a continuous, bounded infinitesimal generator. Thns, we assume 
given an N X N matrix function 



*(z) = (r PQ (z)), 




(24) 


satisfying the conditions 






t pq (z) ^ 0, p * q, t„(z) f£ 0, (p, q = 1, • • • : 


,N), 


(25) 


and 






V 

ErJ = 0, (p = 1, ••• ,N). 

■: 1 




(26) 


Then", for 8z -> +, (/), q = 1, • • • , .V), 






P pl ,(z, z + te) = 1 + t pi ,(z) 8z + o(fiz), 




(27) 


P PQ (z, z + hz) = T» q {z) 8z + 0(&), p ^ q. 




(28) 



If the matrix t is constant then the process M(z) is said to have a 
stationary transition mechanism. 

Stochastic matrix differential equations of the form 

^ = A(I\I(z), z)W(z), W(0) = Y [il/(0)], (29) 
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where W and y are n X m matrices, and A is an n X n matrix, have 
been considered by Morrison. 14 Equations were obtained for calculating 
the stochastic average <W(z)>, in the case that A and y are real, but 
it is easily verified that these equations are still valid for complex 
valued A and y. This may be shown by writing the system (29) in 
real form, and then combining in complex form the equations for 
calculating the stochastic averages of the real and imaginary parts of 
W( 2 ). 

Let Ea, denote the row vector with all N elements equal to 1, and let 

X = (X, , • • ■ , X N ) (30) 

be the row vector of initial probabilities given by (21). Note, from 
(22) and (26), that 

E„X' = 1, £**'(*) = 0, (31) 

where t denotes transpose. Let 

(r kl (z)) p = (r kt (z) | M{z) = a p ) Prob \M(z) = a,,}, (32) 

(fc, I = 0, 1), (p = 1, • • • , N), and introduce the column vectors 

R kl (z) = col ((r kt (*))i , ■■■ , (ru(z)) N ), (/c, t = 0, 1). (33) 

Then, from (9) and (14), 

R kl {z) = (r kt {z)) = E,Rh(2), (/c, t = 0, 1). (34) 

We define the N X N diagonal matrix D(z) by 

D(z) = diag [F(a„ , z)], (35) 

and denote the unit matrix of order N by l N . Then, for the system 
(16) through (19), with N(z) given by (12), it is found that' 4 

^- n = jB(z)(C a R in - CoRo.) + *'(z)Rno , (36) 

^ = jD( 2 )(C ff R,, - C„R„„) + [Ar*I. v + t'(z)]R„, , (37) 

az 

^ = ;D(2)(C' R 1M1 - CRn) + l(Ar + ja)I N + V( 2 )]R 1(1 , (38) 

^ = jvmCRot - CoR,„) + [(AP* + Ar + ja)I, + t'( 2 )]R„ . (39) 
az 

From (20), (21), (30), (32) and (33), the initial conditions are 

R*,(0) = i k i*X, (K I = 0, 1). (40) 
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If the process M(z) has a stationary transition mechanism, and if 
the coupling coefficient (5) depends on z through M (z) only, then the 
differential equations (36) through (39) have constant coefficients. 
We consider a particular example of this in more detail in the next 
section. 

III. RANDOM TELEGRAPH COUPLING 

We now consider the particular case when N(z) = T(z), the random 
telegraph process, 15 which is one of the simplest finite state Markov 
chains. It is an ensemble of square wave functions {T(z)\, such that 
each sample function T(z) can assume only the values ±1. For fixed z, 
a sample function chosen at random will equal +1 or —1 with equal 
probability. The probability p(n, z) of a given sample function changing 
sign n times in an interval of length z is given by the Poisson process 



p(n, z) = *-f- e , 



n\ 



(n- 0,1,2, •••), 



(41) 



where b is the average number of changes per unit length. The process 
T(z) has zero mean and is wide sense stationary, 

(T(z)) = 0, <r(!)T(f)>-«p{-2&|.-r|). (42) 

The states of the process T(z) are a, = 1 and a 2 = —1, and the 
vector of initial probabilities is 



X — (2) 2) — 2^2 — 2"E. 



(43) 



The process has a stationary transition mechanism, with infinitesimal 
generator 



-6 6 

. b -b_ 
Since N(z) = T(z), we have, from (12), 

F(a,,z) = a, = (-I)"" 1 
Hence, from (35), 

"l 



(P = 1, 2). 



D = 



-1 



(44) 



(45) 



(46) 



Note that 



D 2 = I 2 = I, E-c' = 0, ED* 1 = -26ED. 



(47) 
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Let 

Q kt (z) = EDR«(«), (k, I = 0, 1). (48) 

Then, multiplying equations (36) and (39) by E, and equations (37) 
and (38) by ED, we obtain, from (34), (47), and (48), 

^f = KCQ 10 - CoQoi), (49) 

^N = J(C,Q 0l - C Q l0 ) + (Ar* + Ar + jtr)B li , (50) 

dz 



and 



dQot = j(CMu ~ C R Q o) + (Ar* - 2b)Q ol , (51) 

= j(CMoo ~ CoBu) + (Ar + jcr - 26)Q in . (52) 



dz 

dQxa 



dz 

Similarly, multiplying equations (37) and (38) by E, and equations 
(36) and (39) by ED, we obtain 

^f = KCQrr ~ C Q Qoo) + Ar*«o. , (53) 

^§-° = KCjQoo ~ CoQi.) + (Ar + j ff )R l0 , (54) 



and 



dQm = j(CJt 10 - C oJ R 01 ) - 26Q 00 , (55) 

= KCJtoi ~ Coflio) + (Ar* + Ar + j<r - 2b)Q n . (56) 



dQ n 



dz 

From (34), (40), (43), (46) and (48), the initial conditions are 

fi«(0) = i k i* t , QUO) = 0, (k, I = 0, 1). (57) 

Equations (49) through (56), subject to (57), are solved by means 
of Laplace transforms in Appendix A. The Laplace transforms /2 o(*) 
and #„(«) of R 00 (z) and R n (z) are given by (128) and (129), where 
A(s) is given by (130), and d, , d 2 and 77 are defined in (125). Note 
that di and d 2 are linear in the Laplace transform parameter s. Thus 
A(s) is a quartic in s, and we denote the roots of A(.s) = by s k = 
s k (A0, a), k = 1, 2, 3, 4. Inversion of the Laplace transforms in (128) 
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and (129) leads to 

•I 
R 00 (z) = E A ok (A(3, a) exp [s k (A0, a)z] (58) 

and 

4 

Rn{z) = Z A lt (A/3, <r) exp [s t (A/3, <r)z], (59) 

the coefficients A ok and A lk being the residues at s = s k of the expressions 
in (128) and (129) respectively. 

Similarly, the Laplace transforms fi 01 (s) and Ri (s) of R i(z) and 
Ri„(z) are given by (138) and (139), where E(s) is given by (140), 
and e x and e 2 are denned in (135). Note that e t and e 2 are linear in s, 
so that E(s) is a quartic in s. We denote the roots of E(s) = by y. k = 
/i A (A/3, a), k = 1, 2, 3, 4. Inversion of the Laplace transforms in (138) 
and (139) leads to 

•1 
RM = Z Bo*(A/3, <x) exp [ M *(A/3, a)z\, (60) 

and 

4 

Ru>(z) = E B lfc (A/3, CT ) exp [ M *(A/3, «r)a], (61) 

the coefficients B ok and B 14 being the residues at s = n k of the expressions 
in (138) and (139) respectively. 

IV. WHITE NOISE COUPLING 

In this section we investigate the case of white noise coupling, cor- 
responding to (5) and (10), by means of the Ito calculus. We begin 
by considering the vector equation 

j z = Au + Bu£( 2 ), (62) 

where A and B are constant matrices, and £(z) is Gaussian white noise, 
with zero mean and spectral density 1. According to Wonham, 17 the 
Ito differential equation corresponding to (62) is 

d\x = (A + ±B 2 )udz + Budiv, (63) 

where w is a Wiener process. The result is stated for real equations, 
but it may be verified that it holds for complex A, B and u, by con- 
sidering the equations for the real and imaginary parts. 
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The integral equation corresponding to (63) is 

u(z) = u(0) + (A + *B 2 ) f u(f) d{ + B ( u(r) dw({). (64) 

Jo ■» o 

The second integral in (64) is an Ito integral, and it has the property 
that its stochastic average is zero. 18 Thus, 

<u(*)> = <u(0)> + (A + |B 2 ) f <u(f)> dr, (65) 

Jo 

or, in differential form, 

!<u( 2 )> = (A + iB 2 )<u( 2 )>. (66) 

We will also make use of the Ito calculus for the product of differentials, 18 
namely 

(dw) 2 = dz, dwdz = 0, (dz) 2 = 0. (67) 

We now consider the line equations for white noise coupling, so 
that, from (5) and (10), 

c(z) = CVD~o&). (68) 

Hence, from (1), (62), and (66), 



and 



| </„(*)> - -(T, + \C 2 D )(I (z)), (69) 



| (IM) = -(r, + \CD*){m). (70) 

1.9 



These equations agree with those obtained by Rowe and Young. 
From (2), (69), and (70), 

(/,(*)) = U exp [-(r, + hC 2 D )z), (t = 0, 1). (71) 

Also, from (5), (7), (62), (63), and (68), 

dg = -\C 2 D Q g dz + jCVWo 0i dw, ~ 2 * 

dg l = (AT - \C 2 D )g x dz + JCVd'o g dw. 

But, from (14), 

dr kl (z) = dg k (z, A/3 + <r)gf(z, A/3) + g k (z, A/3 + a) dg1(z, A/3) 

+ dg k (z, A/3 + c) dg1(z, A/3). (73) 
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Hence, from (4), (14), (15), (72), and (73), making use of the relation- 
ships in (67), 

dr 00 = D Q [C Csn - \(Cl + Cl)r 00 ] dz 

+ jVK (C a r 10 - O,,,) dw, (74) 

dr 0l = [AT* - \(Cl + Cl)D }r m dz + C C a D r in dz 

+ j\/X^ (C tt r n — C r bo) dw, (75) 

dr l0 = [(Ar + ja) - KCo + Cl)D ]r 10 dz + C C,D r 01 dz 

+ jVWo (C a r 00 - Cor,,) dw, (76) 

dr u = [(AT* + Ar + ja) - \{Cl + Cl)D a \r u dz + C C a D Q r 0Q dz 

+ ; VaJ (C ff r 01 - C n r,o) dw. (77) 

From (9), (14), (63), and (66), it follows that 

^f - -KCo + ODJloo + C C.D R 11 , (78) 

^ - CoCDJloo + [(Ar* + Ar + fr) - *(C5 + C*)Z>o]#n , (79) 

and 

^ = [Ar* - \{Cl + Cl)D ]R 0i + C C a D R 10 (80) 

^ = C„C„A,fl 01 + [(Ar + p) - h{Cl + Cl)D ]R 10 . (81) 

These equations for R kt are consistent with those obtained by Rowe 
and Young. 1,9 The initial conditions are given by (57), and the solutions 
of (78) through (81) are easy to write down. The characteristic roots 
corresponding to (78) and (79) are, from (4), 

X - (Aa + hja) - \{Cl + Cl)D ± [ClClDl + (Aa + \jaf}\ (82) 

and those corresponding to (80) and (81) are 

X = (Aa + li«r) - i(Cn + C;)A, ± [ClClDl - (A/3 + §<r) 2 ]*. (83) 

We now return to the random telegraph coupling and consider a 
case where the coupling is large but has a short correlation length. 
We define the correlation length of the random telegraph process to be 

l. = £ ■ (84) 
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which is just the separation at which the covariance drops to 1/e of 
its peak value. This also agrees with a more general definition of the 
correlation length of a wide sense stationary process, due to C. L. 
Mallows: 10 



t e = lim {'2-k f | S(w) I 2 du \ I S(u) rfw 

o-co l J-n 'I J -« 



(85) 



(88) 
(89) 



where S(o>) is the spectral density function. We set 

C = VbS , C. = VbS„ . (86) 

If we substitute (86) into the expression (130) for A (a), and (140) 
for E(s), as b — > <*>, A(s) has the roots 

8 = (Aa + W - K% + &) ± [(Aa + §j<r) 2 + 8M* + o(|) , (87) 
s = -2b + (Aa + hfr) + H& + 8.) 

± [&S. - (A/3 + i^) 2 ]* + o(|) , 
and #(s) has the roots 

M = (Aa + }jo) - |(& + S.) ± [S S, - (A/3 + |a) 2 ] J + o(|) , 
M = -26 + (Aa + |i«r) + i(& + *S„) 

± [S S, + (Aa + *j<r) 2 ]» + o(±)- (90) 

With the aid of (87) through (90), it is straightforward to calculate 
lim 6 _„ R kt (z), (*, ( = 0, 1). These limits are just the corresponding 
expressions for white noise coupling, if in the expressions for white 
noise we set 

C ( 2 ,D n = So , CD = S a . (91) 

But (91) is consistent with (86), in view of (10), since 

f (T(z + f)T(D) dz = \, (92) 

from (42). 

VI. PERTURBATION THEORY 

We now assume the random coupling is described by a zero mean, 
wide sense stationary process N(z), with an autocorrelation function 
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satisfying (102), but which is otherwise arbitrary. We assume in addition 
that there is a characteristic length t, e.g., the correlation length t c , 
and an e, < e « 1, such that 

CJ = 0(e), CJ = 0(e), Aat = 0(e 2 ), at = 0(e 2 ), (93) 

and in the first instance 



A/3£ = 0(1). 



(94) 



This can best be described as a case with weak coupling, weak atten- 
uation, and narrow fractional bandwidth. 

We consider the asymptotic calculation of R kt (z) for small e. The 
asymptotic results are valid for z/f. = 0(l/e 2 ). We begin by stating, 
in a form which will suffice for our present purposes, some results of 
Papanicolaou and Keller 12 on stochastic differential equations. 

Thus, consider the linear vector stochastic differential equation 



dw 
dz 



= [eA,(z) + e a A 2 ]w(z), w(0) = w , 



(95) 



where Ai(z) is a randon matrix with zero mean, 

<A,(«)> = 0, (96) 

A 2 is a constant nonstochastic matrix, and w„ is a nonstochastic vector. 
Define 

AX = Urn [^ ^ J <A,(2)A 1 (f)> # cfej ■ (97) 

Then asymptotically, for < e « 1, 
d 



j z (w(z)> « e 2 (A,A, + A 2 )<w(z)>, <w(0)> = w„ . (98) 

We now apply the above result to equations (16) through (19) for 
r k ,(z), subject to the assumptions (93) and (94), where Ar is given by 
(4). Thus, we let 

w = (r 00 ,r n , c'^V,,! , e"'^V in )', (99) 

where i denotes transpose. Then, 

-CoT' 4 * CJ*"' 



€A,(2) = jN(z) 



C/ ,w ' -CJ** 






(100) 
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and 



€ 2 A 2 = 





(2Aa + jff) 

Aa 



(101) 



L0 (Aa + ja). 

Note that A, (z) as g iven by (100) satisfies (96), by virtue of (13). 
The calculation of AA , defined by (97), is carried out in Appendix B, 
under the assumption that 



lim [| f' f | P (f) | df J = 0, 
where p(f) is given by (13). Let 

S(«) = r e''" r p(f) dr. 

Jo 

Since, from (9) and (14), R kt (z) = (r kt (z)), it follows from (98), (99), 
(101), (148), and (149) that, asymptotically, 



(102) 



(103) 



dR 0l 
dz 



dR\i 
dz 



-[C 2 S(-A® + ClS(A(3)]R 00 + C C,[S(A/3) + £(-A/3)]fl n , 

(104) 



C o C.[S(A0) + S(-A/S)]ft 



+ [(2Aa + ja) - C 2 S(A(3) - ClS(-A0)]Ru 
and, from (150), that 

^ » [(Aa - jA/3) - (C! + Cl)S(Afi\R ol , 

dz 

^ & [(Aa + jA|S + j<r) - (C? + C$<S(-A0)]B 1O 

The initial conditions are given by (57) . 
Now, 

f cos K) 5(f) df = | / cos («r) 5(f) df = |, 
f°sin(a>f) 5(f) df = 0. 



(105) 

(106) 
(107) 

(108) 
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Hence, from (103), for p(f) = D 5(f) we have S(<a) m 1/2(D ). Note 
that the asymptotic equations (104) and (105) for R 00 and R u are, 
in fact, exact for white noise, as is seen from (4), (78) and (79). However, 
from (57), (106) and (107), for S(u) = 1/2(Z) ), we have asymptotically 

R 0l (z) » totf exp { [Ar* - KCo + Cl)D ]z\ , (109) 

Bio(«) ~ Mo* exp j [(Ar + jo) - |(C" + C:)D ] 2 } . (110) 

It is readily verified that these results are asymptotically consistent 
for ^ z/l ^ 0(l/c 2 ), to lowest order in e, with the exact equations 
(80) and (81) for white noise, under the assumptions in (93) and (94). 
The characteristic roots corresponding to (80) and (81) are, from (83), 



X = ±j(A/3 + \o) 4- (Aa + |;V) - J(Cj + C;)D 4- 



0(e 4 ) 



(HI) 



The above results are valid under the assumptions (93) and (94). 
Suppose now that instead of (94) we have 



(M3)( = 0(e 2 ), 



(112) 



which we refer to as the resonance case. Note if we set b = b /e 2 and 
I = l c = 1/26, the correlation length, the random telegraph case 
considered at the end of Section IV is a special example of this resonance 
case. The resonance case can thus correspond to large coupling, moderate 
attenuation and fractional bandwidth, and small correlation length. 
Corresponding to (16) through (19), we take 



so that, from (93), (95), and (112), we now have 



eA,(z) = jN(z) 



and 










-Co 

c. 








C. 

-Co 



e"A ? = 



(Ar* + Ar + jo) 



o 




-Co 

c, 








Ar* 





-Co 



. 





(Ar + ja)j 



(113) 



(114) 



(115) 
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Note that A x (g) as g iven by (114) satisfies (96), by virtue of (13). 
The calculation of A, A, , defined by (97), is carried out in Appendix B. 
It follows from (98), (113), (115), (152), and (153), with R k( (z) = 
(r kl (z)), that asymptotically 



dR , 

dz 



-(Co + Cl)S(0)R oo + 2C C a S(0)Rn , (H6) 



^ « 2C C a S(0)Roo + [(Ar* + Ar + ja) - (Co 2 + C*)S(0)]fl n , (117) 
dz 



and 
dR i 



dz ~ 



[Ar* - (d + c 2 )S(0)]/e i + 2C c,s(o)# 10 , (lis) 



^ ~ 2C C.S(0)Roi + [(Ar + j<r) - (C{ + Cl)S(0)]R iQ . (119) 

az 

Thus, in the resonance case, the asymptotic equations (116) through 
(119) for R kt (z), (k, I = 0, 1), are exact for white noise, for which 
8(0) = 1/2 (D ), as is seen from (78) through (81). 

APPENDIX A 

We here solve equations (49) through (56), subject to the initial 
conditions (57). The Laplace transform of F(z) is 

F( s ) = [ X e "'F(z) dz. (120) 

Jo 

Taking the Laplace transforms of equations (49) through (52), and 
using (57), one finds that 

s£„ + j'Co&i - jC.Qn = i i% (121) 

(s - Ar* - Ar - j*)& u - jC a Q ol + AOio = *i*t, (122) 



and 



(s+2b - Ar*)Qo, - jtfAi - j'Co^oo , (123) 

(s + 26 - Ar - ja)Q l0 = jCMoo - jCjtxx . (124) 



Let 



d, = s + 26 - Ar*, d 2 = s + 26 - Ar - jV, 

77 = 2Aa + jff. (125) 

Then, substituting for Q 0l and Q i0 from (123) and (124) into (121) 
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and (122), it follows that 

( s + 1 + SK - c "4k + iK = ** - a*» 

Hence 

#00 = {[(« - t?) d, rf 2 + C? d, + Cl d 2 ]i i% + CoCC^ + d 2 )wf}/A(s), 

(128) 
An = [<7oC.(d, + d 2 )io# + (s d, d 2 + Cld 2 + Cl d,)t,tf }/A(«), (129) 

where 

A(s) s [s(fi - n) d, d 2 + (CJ + Cl)s(d 1 + d 2 ) 

- v(Cl d, + C* d*) + (C; - C 2 .) 2 ]. (130) 

Similarly, taking the Laplace transforms of equations (53) through 
(56), and using (57), 

(s - Ar*)£ 01 + jC„Qoo - fCAxx = Mf , (131) 

(s - AT - jo)ft 10 - jC a Qo + jCoQu = iii?! , (132) 

and 

(s + 2b)Q 00 = iC,A,o - j'CoAo, , (133) 

(s+2b- yj)Q u = jCMoi - jC&o , (134) 

where r\ is defined in (125). Let 

e, = s + 26, e 2 = s + 26 - v . (135) 

Then, substituting for Q 00 and Q u from (133) and (134) into (131) 
and (132), 

(a - Ar* + S + ^ oi _ CoCff (l + I^ io m ioif t (136) 
"^'(^ + *)*•« + (• - Ar - *r + ^ + ^)/e 10 = « . (is?) 

From (136) and (137) it follows that 

#,„ = {[(« - Ar - jVte.e, + C 2 e, + Cfcjtotf 

+ C CM + e 2 )idZ\/E{s), (138) 
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and 



£io = \C C a {e, + e 2 )i if + [(a - Ar*Ke 2 + Cfc 2 + CkJMf }/#(•), 

(139) 



where 



£(s) = (s - Ar - ;V)(& - Ar*)e,e 2 + (CS + C 2 ff )s(e, + e 2 ) 

- [Cl&T* + C 2 (Ar + *r)fe - [C 2 (Ar + jV) + C 2 Ar*]e 2 
+ (C 2 - C 2 ) 2 . (140) 



APPENDIX B 



We here carry out the calculation of A, A, , defined by (97). Cor- 
responding to (100) we have, from (13), 

"f o" 



eW^Aitt)) - -p(*- f) 



G 



where 



F = 

and 



[C5r' A ' c '~ n + cy A ' c '" n ] -2C c„ cos a/3(« - r) 
-2C c, cos a/j(* - f) [W^ + c!r' MC -°l 



(141) 



, (142) 



G = 



(143) 



"(C! + Cy*"-" -2C o ae'' A " (2+r) 

L-2C C^- fA/,Cs+n (Co + Of' 4 ' (, ~ r) J 

Now, 

Jo Jo 

= r r e /Bf p(f) <x dz - r (^ - iv*^ #■ (^^ 

Jo Jo *'o 

Hence, for real w, from (102) and (103), 

lim [± f f e iai " S) p(z - f) dt dz] = S(«). (145) 

Z-.W L^ Jo Jo - 1 

Also, 

f f e iu(!+n p(z - f) df d* = / o Z e 2iu ' JT" e-' uf p(f) df & 

= i r z ,«„* r ,-t _ ^"!)p<0 dr. (146) 

2jw J 
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Hence, for real u 5* 0, 

lim [| f* J' e'"<' +f> p( 2 - fl d{ <fe] = 0. 

Tims, from (97), (141) through (143), (145), and (147), 

e'AX = \ Xm ° " 
L Y(A/3). 

where 



(AS)/ = 0(1), 



X(A/3) = 

and 



- [C*S(- A/3) 4- ClS(M)] C C a [S(AP) 4- S(- A/3)] 
L C C ff [S(A/3) + S(- A/3)] - [0(A/3) + Cj8f(- A/3)]J 



Y(A/3) = (C2 4- Cl) 



S(Afl 

. S(-A/3)J 
Now consider A,(z) as given by (114). Then, from (13), 

H 



e"(A,»Ai(f)> = -p(2 - f) 



H 



H = 



where 

"(Cl 4- C;) -2C C, 
. -2C C\ (Cl 4- Cl)_ 

Hence, from (97), (145), (151), and (152), 

H 



e'A.A, = -5(0) 



H 



(A/3K = 0(0 
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(147) 



(148) 



(149) 



(150) 



(151) 



(152) 



(153) 



Note, from (148) and (150), the nonuniform behavior of e 2 A,A! in 
the neighborhood of resonance, which arises from the discontinuity 
of the limit in (147) at co = 0. 
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